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ABSTRACT 

Using high resolution near-infrared spectroscopy with the Keck telescope, we 
have detected the radial velocity signatures of the cool secondary components in 
four optically identified pre— main-sequence, single-lined spectroscopic binaries. 
All are weak-lined T Tauri stars with well-defined center of mass velocities. The 
mass ratio for one young binary, NTTS 160905-1859, is M2/M1 = 0.18±0.01, 
the smallest yet measured dynamically for a pre— main-sequence spectroscopic 
binary. These new results demonstrate the power of infrared spectroscopy for 
the dynamical identification of cool secondaries. Visible light spectroscopy, to 
date, has not revealed any pre— main-sequence secondary stars with masses <0.5 
Mq, while two of the young systems reported here are in that range. We compare 
our targets with a compilation of the published young double-lined spectroscopic 
binaries and discuss our unique contribution to this sample. 

Subject headings: binaries: spectroscopic — stars: pre— main-sequence 



1. Introduction 

The most important characteristic of a star is its mass. Although there have been recent 
advances in determining the masses of 1—2 Mq pre— main-sequence (PMS) stars, there exist 
only a few measurements of the masses of sub-solar mass PMS stars (Simon 2000). As a 
result, theoretical PMS evolutionary tracks have not been well calibrated at the low-mass 
end where uncertainties in the models give rise to scatter as large as a factor of 3 in mass and 
10 in age (Simon et al. 2000). Dynamical measurements of the masses of close binary systems 
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comprised of coeval components with a small mass ratio provide an initial step towards tests 
of the PMS tracks over a broad parameter range. 

In the last decade, much progress has been made towards determining the masses of 
PMS stars. Ghez et al. (1995), Simon et al. (1996), and Steffen et al. (2001) take advan- 
tage of near-infrared (IR) speckle interferometry and the fine guidance sensors available on 

the HST to map the orbits of young binary systems. Using optical wavelength (~5000 A) 
spectroscopy, Mathieu et al. (1989) and Mathicu (1994) identified twelve young, low-mass, 
double-lined spectroscopic binaries (SB2s) and thirteen single-lined systems (SBls). Follow- 
ing this seminal work, other groups have recently detected new PMS spectroscopic binary 
systems. We discuss the sample of published SB2s in §5. 

Our immediate goal is to measure mass ratios and to enlarge the sample to study the 
mass ratio distribution of young stars formed in diverse environments. Precise measurements 
of the mass ratio distribution among young star spectroscopic systems provide critical input 
for theories of star formation. For example. Bate (2000) shows how the density profiles 
of progenitor molecular cloud cores determine the distribution of mass ratios in the close 
(< 1 AU) binary population which forms there. Absolute masses can be obtained from 
SB2s which are eclipsing or from SB2s with known orbital elements determined from high 
resolution mapping of the binary orbits (e.g., Boden et al. 1999). 

SB2s identified by optical wavelength observations tend to have mass ratios close to one. 
As first described in Prato (1998) and further developed in Mazeh et al. (2001), our approach 
to detecting low mass ratio systems is based on the strategy that, at IR wavelengths, the 

flux of the secondary star is a greater fraction of the total flux than at optical wavelengths. 
For example, for two blackbodies of temperatures 5000 and 3000 K, the ratio of the Planck 
functions at 5000 A is about 50, whereas at 1.555 /xm, the ratio is about 4. Mazeh et 
al. (2001) demonstrate the efficacy of this technique, applied to three main-sequence SBls. 
Mazeh et al. (1997) used an intermediate approach, detecting main-sequence binaries with 
mass ratios of 0.57 ± 0.02 and 0.48 ± 0.03 in the R band. With the advent of powerful cross- 
dispersed echelle spectrometers, such as NIRSPEC on the Keck telescope, the IR approach 
has become even more efficient. 

In this paper we describe new observations of four PMS SBls for which we have detected 
the low mass companion by cross-correlation against an extensive new library of spectral type 
standard star templates. The observations and data reduction methods are described in §2, 
the template spectra are presented in §3 and the derived dynamical mass ratios are discussed 
in §4. Section 5 provides a discussion of the results, including a comparison with the complete 
sample of previously published SB2s and a comparison of our mass ratios with the predictions 
of evolutionary tracks in the H— R diagram. A summary appears in §6. Previous reports 
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on applications of this technique have appeared in Prato (1998), Steffen et al. (2001), and 
Mazeh et al. (2001). 

2. Observations and Data Reduction 
2.1. Observations 

2.1.1. NIRSPEC Observations 

Spectra of the target objects were obtained in 2000 June and July and in 2001 January 
and June; specific dates and target properties appear in Table 1. We used NIRSPEC, the W. 
M. Keck Observatory's cross-dispersed near-IR cryogenic spectrometer (McLean et al. 1998, 
2000). NIRSPEC's spectroscopic detector is a 1024 x 1024 ALADDIN III InSb array that is 
sensitive from 0.95—5.5 jim.. In high resolution mode, NIRSPEC yielded a measured spectral 
resolution of R = 24,000 in 9 orders across the H band, centered at A = 1.555yum in order 49. 
Experience has shown that this spectral region contains lines useful for the characterization 
of stars ranging from spectral types of early G through late M, and is reasonably free from 
terrestrial absorption lines (Prato 1998; Mazeh et al. 2001). 

The plate scale is 0."144/pixel in the dispersion direction, and 0."193/pixel in the cross- 
dispersion direction, and the two pixel slit width was 0."288. A 15" nod along the slit enabled 
us to obtain sequential spectra at two positions on the detector for background subtraction. 
NIRSPEC's slit-viewing camera allowed simultaneous monitoring and correction of the ob- 
ject's position on the slit during the spectroscopic exposure. Integration times for a single 
exposure varied from 60 to 600 s, depending on the source brightness. Multiple flat and dark 
frames were median filtered and the results differenced and normalized to create a final fiat 
field. 

Internal Ne, Ar, Xe and Kr arc lamps usually provided wavelength scale calibration and 
dispersion solutions. OH night sky lines, with identifications from Rousselot et al. (2000), 
were also used. A comparison of the night sky line and arc lamp line calibrations for the 
2000 July data showed no difference, but the differences in the 2001 June data between the 
two calibrations was as much as 5 pixels. This drift of the grating position over time on 
that occasion was attributable to vibration caused by the image rotator mechanism inside of 
NIRSPEC. Observations in the adaptive optics (AO) mode (§2.1.2) were never affected by 
grating position drifts since the AO system provides external compensation for field rotation 
and NIRSPEC's internal rotator is not used. 
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2.1.2. NIRSPEC Adaptive Optics Observations 

Because NTTS 155913—2233 is a hierarchical triple, with a tertiary component located 
at a 0."29 separation (Ghez et al. 1993), the 2000 June observations were taken with NIR- 
SPEC behind the Keck II facility AO system (Wizinowich et al. 2000). Magnification by 
the AO reimaging optics alters the plate scale, resulting in pixels a factor of 10.6 times 
smaller. Correspondingly, the two pixel slit width becomes 0."027. The measured spec- 
tral resolution behind the AO system is higher by about 25 %, yielding R — 30, 000. The 
NTTS 155913—2233 tertiary was easily resolved with no overlap of the components' FWHM 
= 0."06 point spread functions. High resolution images, taken immediately before or after 
the spectroscopic integrations, indicated that the Strehl ratio was typically 0.25 in the H 
band. In AO mode, the slow input beam and plane parallel window produce a Fabry-Perot 
type fringing pattern which contaminates the observed spectra. A notch or Hanning filter 
successfully removed most of this artifact during data reduction. 

2.2. Data Reduction 

Two approaches were used for the data reduction. At Stony Brook, we extracted the 
spectra using software originally written by Prato (1998) and adapted to the format of the 
NIRSPEC spectra. At UCLA, spectra were extracted using software developed by S. S. 
Kim, L. Prato, and I. McLean. Comparison of extractions of the same object by the two 
different approaches yielded spectra that were indistiguishable; the locations of the spectral 
fines were identical to within ~ 1 km s~^. Order 49 had no detectable terrestrial absorption 
fines, as expected from prior experience and inspection of the infrared solar atlas (Livingston 
&; Wallace 1991) and order 50 displayed only a few weak terrestrial lines. Unfortunately, 
orders 51, 52, and 53 contained so many strong absorption lines as to be unusable. The 
atmospheric contamination in orders 45—48 was intermediate between these extremes. We 
used only the spectra in order 49, and confirmation from the spectra in order 50, for the 
results reported in this paper. 



3. Template Spectra 

To search for the presence of the secondary in the target spectra, an extensive new suite 
of template standard star spectra were required. Table 2 lists the 23 main-sequence stars 
observed for this purpose. The spectral type appears in column (2), as provided by the Centre 
de Donnees astronomiques de Strasbourg (SIMBAD). Column (3) fists the Jufian Date of 
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the observation and column (4) the derived radial velocity (Mazeh et al. 2001). All but 8 

of the template stars have radial velocities measured by the Geneva (Duquennoy & Mayor 
1991) or CfA groups (Latham 1985). The average difference between the radial velocities 
we measured and those of Duquennoy & Mayor (1991) and Latham (1985) is —0.2 ±0.8 km 
s~^. We therefore adopt ±1.0 km s~^ as an estimate of the systematic uncertainty between 
the zero-points of the velocities we shall derive for the components of spectroscopic binaries 
we observed and their center-of-mass (7) velocities derived by others. 

The excellent sensitivity, resolution, and spectral range obtained with NIRSPEC are 
demonstrated in Figure 1, where the 23 H band template spectra for order 49 are plotted 
with corrections for heliocentric motion and radial velocity applied. The continuum has been 
flattened. The correlation coefficient resulting from cross-correlation of two template spectra 
against each other provides a measure of the stars' similarity. The template spectra have been 
arranged, in Figure 1 and in Table 2, in order of decreasing mutual correlation, from hottest 
to coldest. HD 283750 and GL 275. 2A appear out of order with respect to their spectral 
types, K2 and M4, respectively, as given in SIMBAD. For HD 283750, the discrepancy is 
only 3 spectral sub-classes, however, GL 275. 2A has been incorrectly classified as an M4. Its 
IR spectrum is obviously that of a K star. 

Based on Poisson statistics, the signal to noise ratio of the spectra shown in Figure 
1 is typically greater than 100; most of the observed features are real. The G through 
early M type standard star spectra are dominated by atomic lines, predominantly Fe I, 
Ni I, Si I, and Ti 1. Beginning around spectral type K5, OH lines appear and increase in 
depth. A pronounced transition is evident between the M3 and M4 subclasses; the earlier 
type standards generally show a clearly definable continuum, whereas objects later than M3, 
dominated by water features, do not display any region clear of lines from which a continuum 
might be definitively characterized. In addition to H2O, late M star spectra in this region 
contain strong OH hnes and the CO Av — 3 bandhead at ~1.558 //m. 



4. Results 

We selected the targets for our study from the sample of SB Is described by Math- 
ieu (1994) and here report on the results for four systems, NTTS 155913—2233, NTTS 
160905—1859, Parenago 1771, and Parenago 1925; observations of other systems are still in 
progress. A summary of the properties of the four target SBs appears in Table 1. Each of 
these binaries was detected as an X-ray source by the Einstein observatory (Ku et al. 1982; 
Walter 1986; Mathieu et al. 1989; Gioia et al. 1990); all are weak-lined T Tauri stars with 
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deep photospheric absorption lines. The 7 velocities have already been well determined, from 
optical observations, for these four systems (Mathieu 1994). Figure 2 shows representative 
order 49 spectra of the targets as observed without correction for heliocentric velocity. The 
continuum has been flattened for each spectrum. 

We analyzed the spectra, following the procedure described in Mazeh et al. (2001), with 
the TODCOR algorithm, which uses a two-dimensional cross-correlation for the analysis 
of composite spectra (Zuckcr & Mazeh 1994). Table 3 presents the results for each target: 
columns (1) and (2) list, respectively, the Julian date and the orbital phase of the observation. 
TODCOR analysis identifies the templates that provide the best fit to the primary and 
secondary, their velocities, vi and V2, given in columns (3) and (4) of Table 3, and an 
estimate of their flux ratio (see Table 4) . The 7 velocity is known for our targets from prior 
analysis as SB Is (Mathieu 1994), so the mass ratio of the components, q — M2/M1, column 
(5) of Table 3, follows from their velocities relative to ^, q — {vi — 7)/(7 — V2) (Mazeh et al. 
2001). Values for 7 appear at the top of each section of Table 3 adjacent to the target names. 
For all 4 of our targets, the velocities listed in Table 3 are the averages of the analyses of 
order 49 using different template stars. Preliminary results from order 50 were consistent 
with these. The errors in the component star velocity determinations were propagated to 
determine a final uncertainty in the mass ratio. 



4.1. NTTS 155913-2233 

This NTTS is a hierarchical triple, consisting of an SB with a period of 2.42 days and 
a tertiary component located at 0."29 separation (Mathieu et al. 1989; Ghez et al. 1993; 
Walter et al. 1994). Our observations of this multiple were obtained using the AO system 
to separate the SB and the tertiary. TODCOR analysis identified GL1094, rotationally 
broadened to 15 km s"-*^, as the best-fitting template for the primary. Its spectral type, K5, 
agrees with the determination of Mathieu et al. (1989) for the SB seen in visible hght (Table 
1). GL436, also rotationally broadened to 15 km s~^, provided the best fit to the secondary, 
with a flux ratio of 0.39 ±0.03 for two of the observations. On UT 2000 June 11, however, 
TODCOR identified a flux ratio of 0.79 ±0.03, although the mass ratio was identical. Since 
NTTS 155913—2233 has a relatively short period and the components are probably only a 
few solar radii apart, the stars potentially fill their Roche lobes (Mathieu et al. 1989), and 
hence are perhaps distorted, producing variation in their fiux ratio with phase. 

The upper left two panels of Figure 3 are cuts through the two-dimensional correlation 
function, results of the 2000 June 6 data TODCOR analysis. The primary and secondary 
correlation peaks, located at ~ ±43 and ~ —74 km s~^, upper and lower panels, are clearly 
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seen. The correlation plot for the secondary also shows a residual bump at ~ +43 km s~^, the 
result of the imperfect match of the primary template to the primary of NTTS 155913—2233. 
Table 3 lists the component velocities measured in our three observations of the system. The 
derived mass ratio, using the 7 velocity measured by Mathieu et al. (1989), is 0.63 ± 0.02. 
This value is the average of two independent observations of the system. The value of q in 
parenthesis has a large uncertainty because it was measured when the component velocities 
were near the 7 velocity. Therefore, we did not include it in the average. 

4.2. NTTS 160905-1859 

The maximum correlation was obtained for NTTS 160905—1859 using the K2 star HD 
283750, rotationally broadened to 25 km s~^, as the primary template. This spectral type is 
consistent with the determination of Mathieu et al. (1989) for the primary. For our highest 
signal to noise spectra, obtained without the AO system, the TODCOR correlation was 
maximized by using GL 402, an M4, as the secondary template, also broadened to 25 km 
s~^. An average flux ratio of ~ 0.24 was derived. The range of flux ratios found with 
different templates and derived from distinct observations was 0.15 — 0.35. Because of the 
lower throughput and hence signal to noise ratio in the AO measurements, it was not possible 
to extract secondary star velocities from these data. The upper right panels of Figure 3 show 
the results of the correlation analysis of the UT 2001 June 2 data. Table 3 summarizes the 
measured velocities. The ratio derived using the 7 velocity from Mathieu et al. (1989) is 
0.18 ± 0.01, the lowest yet measured for a PMS spectroscopic binary. 

4.3. Parenago 1771 

GL1094, rotationally broadened to 15 km s^^, gave the largest correlation as the primary 
of Parenago 1771. Its designated spectral type is K5, in good agreement with the K4 type 
given for the SBl by Mathieu (1994). GL436, broadened to 15 km s~^, with an average flux 
ratio of 0.26 with respect to the primary, gave the best correlation for the secondary. For 
different templates, the range of flux ratios was 0.22 — 0.32. The lower left hand plots in 
Figure 3 show the correlation for the primary (upper panel) and secondary using the order 
49 spectra. In the lower panel, the correlation peaks for the secondary at ~ 64 km s"^ and 
for the residual of the primary at ~ 14 km are nearly equal in amplitude, the result of 
an imperfect match between the main-sequence primary template and the primary in this 
PMS binary. The derived mass ratio is 0.34 ±0.03. It is necessary to conflrm this mass ratio 
because it is based on only one observation. 
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4.4. Pcirenago 1925 

TODCOR analysis of the Parenago 1925 system identified HD 283750 and GL669B, 
each rotationally broadened to 10 km s~^ and with an average flux ratio ~ 0.23, as the 
best fitting primary and secondary. For different templates, the range of flux ratios was 
0.12 — 0.35. The assigned spectral type of HD 283750 is K2, in good agreement with the 
classification by Mathieu (1994) of the primary as a K3. The lowest right panel of Figure 
3 clearly shows the correlation maximum associated with the secondary, but also shows a 
comparable residual peak associated with the primary. The mass ratio in this system is 
0.31 ± 0.04. As for Parenago 1771, since this result is based on only one observation it must 
be confirmed. 



5. Discussion 

To illustrate the sensitivity of the IR technique. Figure 4 combines our sample with 
the IR detected SB2 NTTS 045251+3016 (Steffen et al. 2001) and 26 optical wavelength 
detected SB2s. From Figure 4 it is evident that the small sample presented in this paper 
includes the lowest mass ratios ever measured for PMS SB2s. Table 4 lists the entire sample 
of PMS SB2s plotted in Figure 4, including the location of the binaries, column (2), their 
mass ratios, periods, flux ratios, and primary masses (see below), in columns (3), (4) (5), 
and (6), and the references used for these data, column (7). 

Figure 4 displays the well-known bias of SB2s detected in visible light toward mass 
ratios of one. Our results start to fill in the low mass ratio part of the observed distribution. 
A concentrated effort to observe a complete sample of SBs will be required to determine 
the true distribution of PMS SB2 mass ratios. It is interesting to notice (Table 4) that the 
five SB2s contributed so far by the IR technique span the same range of periods identified 
among the optically detected SB2s. A complete sample of SB2s will therefore enable tests 
of theoretical predictions of the mass ratio distribution as a function of binary separation 
(e.g.. Bate and Bonnell 1997). 

High angular resolution imaging studies of young binaries reveal a broad range of 
component mass ratios, down to the regime of brown dwarfs, as inferred photometrically by 

comparing the location of the binary components in the H— R diagram with PMS evolution- 
ary tracks (White et al. 1999; White & Ghez 2001; Woitas et al. 2001). However, differences 
among the theoretical calculations of PMS evolution produce considerable scatter among the 
inferred mass ratios (e.g.. Figure 7 of Woitas et al. 2001). Dynamical measurements offer the 
means to determine the mass ratios and their distribution precisely and thus to test specific 
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predictions of calculations of binary formation (e.g., Bate 2000). 

Except for the eclipsing systems, indicated in Table 4, and the spatially resolved 
NTTS 0455251+3016 (Steffen et al. 2001), we cannot derive dynamical values for the com- 
ponent masses of the SB2s in Figure 4 because their orbital inclinations are not yet known. 
In order to compare our IR SB2 secondary masses with their minimum values set by the 
mass function, and to learn the mass range currently sampled by dynamical techniques, we 
estimate the primary masses, and thus M2 using the dynamical mass ratios (Table 5). The 
procedure we use to obtain Mi is described in detail in Appendix A. The Mi values for the 
entire SB2 sample are included in Table 4, column (6). 

For the 4 IR SB2s described in this paper, columns (3), (4), and (7) of Table 5 
provide estimates for Mi, M2, and the minimum mass, M2^-„, derived by setting the orbital 
inclination to 7r/2 in the mass function f(M) of Mathieu (1994), column (6). The uncertainties 
in Mo and M2 include the uncertainties in the mass ratio and in our estimate of Mi. In 
all cases, M2 is consistently greater than its predicted minimum value. 

In Figure 5 we plot mass ratio as a function of estimated Mi. Figure 5 shows that 
smaller secondary masses, as well as smaller mass ratios, are more readily identified with 
IR spectroscopy; all the companions detected by the IR technique have M2 < 1 Mq. In 
general, few primary star masses less than 1 M© and no primary star masses less than 0.7 
Mq have been detected in PMS SB2s. This is illustrated by the gap on the left hand side 
of Figure 5 and is probably the result of a selection effect for brighter systems among SB 
surveys. The total sample shown in Figure 5 is weighted towards PMS stars with little or 
no circumstellar material, although spanning a wide range of ages, from ~ 10^ yr through ~ 
10'^ yr. Young stars which lack circumstellar envelopes and disks do not experience veiling 
of their photospheric absorption lines, and hence provide easier targets for the identification 
of spectroscopic binaries. 

The essential result of TODCOR analysis is the mass ratio of the components. This 
value is robust because it is determined by cross-correlation of the many spectral lines in 
common between the target and templates. The analysis also provides estimates of the 
primary and secondary spectral types and their H band fiux ratio, which give the maximum 
correlation with the target. The templates almost certainly do not match exactly to the 
spectra of the stars in the target binaries because the surface gravities of PMS stars are 
lower than those of main-sequence stars, and because the templates and the PMS targets 
may differ in metallicity. Similarly, H band flux ratios do not necessarily represent accurate 
component luminosity ratios. Nonetheless, it is of interest to investigate the extent to which 
the values of M2 / Mi implied by our dynamical results are consistent with PMS evolutionary 
tracks in the H— R diagram, presented in Figure 6. 
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Our approach to estimating the effective temperature, Te//, and luminosity for each 
object, and the associated uncertainties, is described in Appendix B. We use the tracks of 
Palla & Stabler (1997) because they are representative of PMS evolution, similar to the 
results of Baraffc et al. (1998) and Siess et al. (2000), and span the range of masses of our 
sample objects. Figure 6 shows that, to within 1 c, the secondary masses implied by the mass 
ratios (Table 5) are consistent with the secondary masses estimated from the tracks (column 
(5) of Table 5), except for Parenago 1925. Only for NTTS 160905-1859 do the components 
lie on distinct isochrones; the other systems appear to be coeval to 1 a. It will be possible 
to derive absolute, dynamical masses for these objects within the next decade as ground 
and space based interferometers become available for mapping their orbits. These mass 
determinations will be precise to within a few percent or less, a requirement for meaningful 
testing of models. Reliable luminosity estimates, however, cannot be obtained by dynamical 
means and will require not only spatially resolved observations, but also multi-wavelength 
data in order to obtain precise values. 



6. Summciry 

We sucessfuUy observed 4 PMS SBl systems using IR spectroscopy to find cool com- 
panion stars not detected at optical wavelengths. A suite of spectral type standard stars 
provided templates for cross-correlation analysis and yielded a spectral library of early G 
through late M type stars. We have identified the smallest mass ratio PMS SB2s yet ob- 
served, demonstrating the power of this approach. The binary components we have studied 
cover a large mass range, from 2 M© primaries to 0.2 M© secondaries. 

To provide a comparison to our sample and a useful reference, we compiled as complete 
a list as possible of PMS SB2s and their main properties. We estimate primary masses of all 
known PMS SBs (Table 4) and quahtatively use these values, together with the precise mass 
ratios, in Figure 5 to illustrate the extremely low secondary star masses discovered with our 
IR observations. In addition, we place both components of the PMS SBs systems studied 
in our sample on the H— R diagram to compare the evolutionary tracks with our dynamical 
results. To within 1 cr, most of these results are consistent, however, mass estimates based 
on placing the stars on the H— R diagram are still approximate because the spectral types 
and luminosities of the component stars are incompletely characterized. 

The application of optical wavelength spectrometers on 1—2 m class telescopes for the 
identification of SB Is, to be subsequently examined in the IR, is a practical technique for the 
measurement of large samples of PMS SB2s. Although such a study is magnitude limited, 
and thus does not sample the smallest mass primaries, the detection of SBs [not SB2s) is 
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unbiased with respect to the mass ratio. We anticipate that, for a given star forming region, a 
complete survey of SBs at optical wavelengths with IR follow up for large mass ratio systems 
would provide the true mass ratio distribution for the smallest separation observable systems, 
to the limiting magnitude of the initial optical light survey. Because of the multiple epochs 
of observations required, as well as the limited time available at large telescopes, progress 
on this research is inevitably slow-paced. As we proceed with new confirmed observations, 
we plan to pubhsh these results as they become available. 
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A. Primary Mciss Determinations for Figure 5 

Primary masses for fourteen of the binaries in Table 4 were found in the literature 
(column 7). Among these, only for NTTS 045251+3016 (Steffen et al. 2001) is there a 
published value for the fiux ratio. To obtain the primary masses for the remaining systems, 
we did the following. The total luminosity of the SB2 was obtained from the literature or 
calculated from the total K band magnitude, as described in Simon & Prato (1995). The 
primary luminosity was then determined from the fiux ratio, when available. Otherwise, if 
the mass ratio of an SB2 was within 3% of unity, a component bolometric luminosity ratio 
of 1.0 was assumed. For the remaining three SB2s, we used an approximate fiux ratio based 
on the empirical mass ratio — Am.K relation shown in Figure 5 of White & Ghez (2001). 

For cases in which no component spectral type data was available, we assumed that 
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the spectral type of the unresolved system was equal to the Te// of the primary. For the 4 
systems described in this paper, we estimated Teff from the best fit primary star template 
(§4) and the spectral type to temperature conversion found in Luhman (2000) (see Appendix 
B). 

Primary star masses were derived from the tracks of Palla & Stabler (1999) and 
Baraffe et al. (1998) for the lower mass objects. Objects plotted on both sets of tracks usually 
had comparable masses to within 10 %; for consistency we used Palla & Stabler (1999) for 
all estimates. We estimate that the uncertainties in the primary masses determined from 
the H— R diagram and listed in Table 4 are about 10%. 

B. Te// and Luminosity Estimates for Figure 6 

Spectral types for the 8 individual stars in our SBs sample were determined from 
the spectral types of the templates which gave the maximum correlation. Because of the 
caveats in the accuracy of this approach (sec §5) we assigned an uncertainty to these results 
of two spectral subclasses for secondary stars. In all cases the primary star spectral types 
derived from the TODCOR analysis were the same as previous estimates of the spectral types 
for these stars, or else within one spectral subclass. Therefore, we assigned a one subclass 
uncertainty to the primary star spectral types. Typically the same component velocities were 
determined using a range of template spectral types. To convert to Tg//, we used Figure 5 
of Luhman (2000). 

To derive component luminosities we used the 1.555 pLm flux ratio found by TODCOR 
and apphed this proportionately to the total luminosity. For the PMS stars of the age 
and mass range considered here, their H band flux is approximately proportional to their 
luminosity. We estimate one source of uncertainties from the range of flux ratios produced 
by TODCOR for different templates (see §4). For example, for Parenago 1771, an M3 
template gave the best correlation for the secondary star, yielding a velocity of 64.5 ±1.5 
km s~^ (Table 3). An M7 template yielded only a 1 % difference in velocity, yet a ~25 
% difference in the flux ratio. Thus, although the determination of the velocity is robust, 
the characterization of the flux and hence luminosity ratio is uncertain. Because these error 
estimates are approximate, we use a ~15 % distance uncertainty, or else the uncertainty 
produced by the range in flux ratios, whichever is larger. 

Other sources of error might arise from the proximity of the stars in very short period 
SBs, which may distort their structure. The light ratio variation (§4.1) detected in NTTS 
155913—2233 (period = 2.4 days) may indicate this effect. Such systems are likely to be less 



-13- 



useful for tests of the PMS evolutionary tracks. 
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Template Spectra (NIRSPEC Order 49) 
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Fig. 1. — The NIRSPEC template spectra for order 49. The spectral continuum has been 
flattened for each template. 
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Fig. 2.— Order 49 NIRSPEC spectra of NTTS 155913-2233, NTTS 160905-1859, Parenago 
1771, and Parenago 1925, observed respectively on UT 2000 June 6 (AO), UT 2001 June 
1, and, for the latter two, UT 2001 January 5 (AO). The Parenago 1771 and Parenago 
1925 spectra have been smoothed with a two-pixel boxcar. No heliocentric or radial velocity 
corrections have been applied. The spectral continuum has been flattened for each object. 




Fig. 3. — TODCOR correlation results for the spectra in Figure 2. For each binary, the 
upper panel shows the two-dimensional cross-correlation as a function of the radial velocity 
of the best fitting primary template, holding the secondary velocity fixed. The lower panel 
shows the cross-correlation as a function of the secondary template velocity. 
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Fig. 4. — Histogram of currently known SB2 mass ratios. The systems shown, their mass 
ratios, and other data are hsted in Table 4. 
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1 2 3 4 5 6 7 

Fig. 5. — Mass ratios as a function of primary star mass. Triangles are SB2s discovered at 
optical wavelengths and filled circles are SB2s discovered in the IR. Primary masses were 
taken from the literature or calculated as described in the text, §5. The dot-dash fine shows 
the curve for a constant M2 =1.0 Mq and the dotted line for a constant M2 =0.5 Mq. 
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Fig. 6. — Both components of the target systems plotted on the H— R diagram, with the 
evolutionary tracks of Palla & Stahler (1999) superimposed. The plotting symbols are as- 
signed as follows: open triangles, Parenago 1925; pentagons, Parenago 1771; squares, NTTS 
160905-1859; filled triangles, NTTS 155913-2233. 



Table 1. Summary of Sample Properties and Observations 



Object 


R.A. (2000.0) 

{h, m, s) 


Dec. (2000.0) 


HtotaZ 

(mag)*^ 


Spectral 
Type^ 


Period 
(days)"" 


UT Date and Mode 
of observation 


Parenago 1771 


05 35 09.75 


-05 23 26.8 


9.6 


K4 


149.5 


2001 January 5 (AO) 


Parenago 1925 


05 35 18.34 


-05 22 37.7 


8.6 


K3 


32.94 


2001 January 5 (AO) 


NTTS 155913-2233 


16 02 10.5 


-22 41 29 


8.2 


K5 


2.42378 


2000 June 7, 10, 11 (AO) 


NTTS 160905-1859 


16 11 59.33 


-19 06 52.40 


8.2 


Kl 


10.400 


2000 June 9, 11 (AO) 



2000 July 14 (non-AO) 
2001 June 1, 2 (non-AO) 



Data from Mathieu et al. (1989). 



Table 2. Template Library 
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Object 


Type^ 
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1 Q1 fi Qn 
iy iu.yu 


'^Q fi 
oy .u 
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K5 


1703.03 


-67.1 


BS 8086 


K7 


1707.13 


-65.0 


GL 763 


MO 


1715.11 


-60.5 


GL 15A 


M1.5 


1707.14 


12.2 


GL 436 


M2.5 


1915.14 


10.0 


GL 752A 


M3 


1705.03 


37.4 


GL 213 


M4 


1917.96 


106.6 


GL 402 


M4 


1915.04 


-1.1 


GL 669B 


M4.5 


1705.99 


-33.8 


LHS 292 


M6.5 


1915.06 


1.3 


GL 644C 


M7 


1943.17 


15.0 


LHS 2351 


M7 


1942.97 


-14.4 


LHS 2065 


M9 


1942.94 


-8.4 



""As given in the SIMBAD database. 



Table 3. Mass-Ratio Determination 



JD vi 


V2 




(2450000+) Phase (km s"^) 


(km s~^) 


q 



Parenago 1771 7 = 25.6 ± 0.2 km s"^ 
1914.86 ... 12.5 + 1.0 64.5 + 1.5 0.34 + 0.03 



Parenago 1925 7 = 25.9 + 0.3 km s"^ 
1914.94 ... 34.7 + 1.0 -2.2 + 1.5 0.31 + 0.04 



NTTS 155913-2233 7 = -2.3 + 0.7 km s'^ 
1702.96 0.22 43.4 + 1.1 -74.4 + 1.5 0.63 + 0.02 
1705.93 0.45 -41.8 + 1.1 61.3 + 1.6 0.62 + 0.03 
1706.92 0.86 0.4 + 1.0 -5.8 + 1.1 (0.77 + 0.45)^ 

^=0.63 + 0.02 



NTTS 160905-1859 7 = -6.4 + 0.2 km s'^ 



1704.91 


0.73 


-22.3 ± 1.1 






1706.94 


0.93 


-15.6 + 1.1 






1739.89 


0.10 


8.4 + 1.0 


-93.3 + 2.1 


0.17 + 0.01 


2061.95 


0.07 


7.0 + 1.0 


-79.7 + 3.2 


0.18 + 0.02 


2062.91 


0.16 


10.8 + 1.0 


-99.3 + 2.3 


0.19 + 0.02 



^=0.18 + 0.01 



^Not used in q. See §4.1. 



Table 4. PMS Double-Lined Spectroscopic Binaries 
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0.18 


10.4 
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2 


AK Sco 
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0.99 


13.61 


~ 1 


1.6 


4. 7, 8 


HD 155555 
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0.93 


1.68 


0.8 


1.1 


4, 7, 8, 11 


V4046 Sgr 


isolated 


0.94 


2.42 




0.86 


18 


TY CrA'^ 


Corona Australis 


0.52 


2.89 




3.16 


2 


EK Cep^ 


isolated 


0.55 


4.43 




2.03 


2 



Table 4 — Continued 



Object Location q P (days) Flux Ratio Mi (Mq) Refs^ 



^References are for data used to estimate Mi, not for discovery pa- 
pers. 

'^Flux ratio is average of a range of values determined from use of 
several H band templates in the cross-correlation analysis. See §4. 

'^Eclipsing systems. 

References. — (1) Covino ct al. 2001; (2) Palla & Stabler 2001; (3) 
Lee et al. 1994; (4) Luhman 2000; (5) Mathieu et al. 1997; (6) Steffen 
et al. 2001; (7) Mathieu 1994; (8) PaUa & Stabler 1999; (9) Covino et 
al. 2000; (10) tbis work; (11) Wbite & Gbez 2001; (12) Torres et al. 
2000; (13) Padgett & Stapelfeldt 1994; (14) Torres et al. 1995; (15) 
Prato et al. 2001; (16) Walter et al. 1994; (17) Guentber et al. 2001; 
(18) Quast et al. 2001. 



Table 5. Detailed Sample Properties 















f(M)<i 


M2 . 


Object 




qa 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


Parenago 1771 




0.34±0.03 


1.4±0.5 


0.48±0.18 


0.34±0.17 


0.024 


0.43±0.10 


Parenago 1925 




0.31±0.04 


2.3±0.6 


0.71±0.21 


0.15±0.22 


0.011 


0.44±0.08 


NTTS 155913- 


-2233 


0.63±0.02 


1.0±0.2 


0.63±0.13 


0.34±0.20 


0.064 


0.53±0.06 


NTTS 160905- 


-1859 


0.18±0.01 


1.2±0.2 


0.22±0.04 


0.16±0.21 


0.0042 


0.20±0.02 



^Prom Table 3. 
''Prom Table 4. 

" qxMi 

^Prom Mathieu (1994). 



